Quantum corrections to the thermodynamic potential of interacting Bosons in a trap 
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We calculate the quantum corrections of the thermodynamic quantities of a system of confined 
Bosons at finite temperature. Systematically quantum corrections are written in a series of h, which 
is convergent when kT is much larger than the spacing between energy levels of the system. We apply 
this method to calculate analytically the thermodynamic potential of a weakly interacting Bose-gas 
confined in 3-d harmonic oscillator potential. For large number of particles, quantum corrections 
become small, and contribute to the finite size corrections to scaling. 
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There has been renewed interest in the Bose-Einstein 
condensation(BEC) after its experimental demonstration 
by Anderson et al Q| . It has now become possible to mea- 
sure the relevant thermodynamic quantities of a system 
of weakly interacting Bosons in a magnetic trap ||,|| . In 
these experiments, the temperature dependence of the 
condensate fraction and the release energy are measured. 
Theoretically, the thermodynamic properties of homoge- 
neous Bose gas was studied by several authors starting 
from London The thermodynamic properties of the 
inhomogeneous Bose gas was studied by Legget et al ||D ■ 
After the recent experimental developments, there is a 
new impetus to understand the thermodynamic prop- 
erties of interacting inhomogeneous Bose gas. Within 
the local density approximation(LDA), Stringari et al |7]] 
have calculated the thermodynamic quantities numeri- 
cally and shown scaling properties of strongly interacting 
bosons in a trap and since then several calculations have 
been reported by various authors I^Q- In these papers 
the thermodynamic properties of inhomogeneous gas has 
been calculated within the LDA, neglecting the quantum 
corrections which may be important. 

In this paper we develop a systematic semiclassical H- 
expansion for thermodynamic quantities of an interacting 
bose gas which takes into account the quantum correc- 
tions. Wc follow a method analogous to the extended 
Thomas-Fermi method(ETF) which has been success- 
fully applied in the case of finite fermion systems, like 
nuclei, atoms, clusters etc. It is well known that the 
semiclassical expansion of the smooth part of the quan- 
tum density of states can be expanded in powers of h. 
The leading term in this expansion is the Thomas-Fermi 
level density. Systematic Ti corrections to the leading ap- 
proximation may then be obtained using the well known 
Wigner-Kirkwood(WK) expansion (H}. In the case of fi- 
nite fermion systems this is what is usually referred to as 
the extended Thomas- Fermi approximation [Tl| ]. 

We first develop a systematic Ti expansion using the W- 
K method pOj . Consider the single particle Hamiltonian 
H for a system of particles, which satisfies the equation, 



Hip n = E n ip ri 



(1) 



where E n are eigenvalues and ip n are corresponding wave 
functions. The canonical partition function of the system 



is defined as 



(2) 



where (3 is, for our purposes, a parameter with the di- 
mension of inverse energy. We first briefly review the 
method as applied to the finite fermion system jyj. To 
extract thermodynamic quantites of the finite fermion 
system from the canonical partition function, a weight 
function is defined aswj? = sin(f/3T) ■ ^ using this weight 
function one can write the grand potential of the system 
in the following way, 



kTq(T, z) = fiN -E + TS 



x Z (P)w F 
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In integral representation, 



kTq(T,z) 
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where the contour is closed from left, enclosing all the 
negative poles. Summing all residues we obtain the 
known result |l2j , 



q (T,z) = J2(.l + zz- EJkT ), 



(5) 



where z = e" ,i is the fugacity. We may extend this 
method to Bosons also, by taking a suitable weight func- 
tion, 
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- C0S(7T/?T)7T/3T 

sin(7r/3T) ' 
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The grand potential of an ensemble of bosons can be 
written in terms of the following integral, 



z '" <■ J n—ioo H 



(7) 



rj — too 



The line of integration is choosen in such a way that all 
the positive poles (excluding zero) are to the right side 
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of the line. Since the chemical potential /i < E n for 
bosons, we close the contour from right, and summing 
all the residues at the positive poles inside the contour, 
we obtain the well known result |l2|] , 



q(z 7 T) = -Y,log(l-ze- E ^ kT ). 



(8) 



From the above thermodynamic potential we can calcu- 
late all thermodynamic quantities. 

To do a systematic Ti expansion of Bosons, we start 
from high temperature limit. Neglecting zero point en- 
ergy and the occupation number of ground state, we can 
do W-K expansion of canonical partition function. For 
the semiclassical expansion of the partition function, it is 
most convenient to take plane wave basis for calculation 
of the canonical partition function. 

Z (f3) = -L J d 3 P J d 3 re~ l P- f/h e- pH e^ /h (9) 

To evaluate the partition function, we write, 



(10) 



where H c i is the classical hamiltonian H c i — + V(r) 
and u obeys the Bloch equation, 



du 
dp 



Hu = 



with the boundary condition 

lim u = e l P- F/h . 



(11) 



(12) 



Substituting u in the above equation, we obtain the fol- 
lowing equation for w. 

d ^ = -in[^v )w -^ip.vv) ] + 

h 2 

— [f3 2 (VV) 2 w - (3(y 2 V)w + V 2 w - 2/?(W.Vw;)] (13) 

Ziyi 

The above equation can be solved order by order in K, 
by expanding w as a power series in h. 



w = i + hwi + h w 2 



(14) 



Upto the order h , canonical partition may then be writ- 
ten as fill, 



Z wk {j3) 
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(15) 



where V(r) is the effective single particle potential. In- 
serting this expression for canonical partition function in 



equation (7), we obtain the grand potential of the normal 
state. 
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After doing the p integration we get the following form 
of the density of the grand potential upto order h 2 . 
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F{r) = ^\9 5/2 (z) - —-VV(r)— g 5/2 (z)] (17) 



A 3TO/A ' 2AM y 'd^ 2 
where the g n (x) and A are defined as, 

2nh , , 

= ; g n (x)=^ 

i=l 



A = 



V2vrMfcT 



(18) 



and the effective fugacity is z = ze~@ v ( r \ 

To check the correctness of the above formalism, we 
first apply it to a system of bosons in 3-d isotropic 
harmonic oscillator (h.o) confinement. Following our 
method, grand potential of the normal state is given by, 

kT 1 

q(z,T) = (j-) 3 l94(z) - -{Phu) 2 g 2 {z)\ (19) 

From this grand potential we can derive the number of 
particles and the energy of the normal state. 



N e = z{A 
oz 

kT 1 

= ( J -ri9 3 (z)-~(f3Tiu 1 ) 2 g 1 (z) 



(20) 



kT 1 

= kT( — ) 3 [3g 4 (z) - -((3hu,) 2 g 2 (z)} (21) 

Above quantities may also be derived by using density of 
state of the system. Exact canonical partition function 
of the 3-d h.o is given by, 



Z(P) 



1 



[2sinh(^)]3' 



(22) 



Density of states can be calculated by taking inverse 
laplace transform of the canonical partition function with 
respect to (3. Also the same results may be obtained by 
Eular-Maclaurin summation metod |9) . From the above 
calculation we have seen that, thermodynamic quantities 
may be expanded in a series of a dimensionless parameter 
[3huj. This series will converge when f3Tiuj < 1. However 
at low temperatures, {3htv ~ 1, and Wigner-Kirkwood 
expansion breaks down. At these temperature, only few 
low energy states are occupied, so energy levels can not 
be taken as continuous. 
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In the systems which show macroscopic occupation of 
ground state at finite temperature, chemical potential 
of the system tends to zero. As z — > 1, the higher or- 
der terms in h expansion show infrared divergence. To 
regulate such divergences fugacity of the system can be 
replaced by ze~ AE / kT , where AE is the natural energy 
gap of finite size systems |l3| . This is equivalent to in- 
troducing a infrared momentum cutoff *=s 1/L, where L 
is the system size. 

Having checked the formalism, we now apply above 
method to real physical situation that of a non ideal 
Bose gas. In recent experiment Ensher et al [0 measured 
the thermodynamic quantities of 40000 Rb 8 atoms, for 
which s-wave scattering length is rs lOOao, where ao is 
Bohr radius. 

Within Hartree-Fock theory of bosons, non-condensate 
part satisfies the following Schrodinger's equation ||, 



2M 



V 2 + -Mw 2 r 2 



■2up(r)]ip i (r) = e i ip i (r), (23) 



and the condensate part satisfies Gross-Pitaevskii equa- 
tion, 



ff 1 

[ "2M V2 + 2 + 2u P nc ^ + u Pcir)]<f>(r) 



ea<f>(r), 



(24) 



where p n c{r) and p c {f) are the densities of non- 
condensate and condensate part, and p = p c 
interaction strength is given by, u — l,T " " 



+ p nc . The 
Neglecting 

the kinetic energy of the condensate, the lowest energy 
eigenvalue may be approximated by, 



M 



fo 



^hco(l5N Q j) 2 ^ + 2u Pnc {0). 



(25) 



Further within the same approximation condensate den- 
sity can be written as, 
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Pc (r) = -(e 1 -V(r)), 



(26) 



where e% = ^huj(15N f) 2 / 5 . The above approximation 
is valid for large number of condensate atoms Nq and for 
strong repulsive interaction. This approximation breaks 
down within a small region near the critical tempera- 
ture, where the number of condensate atoms become very 
small. For the non-condensate part we may use the for- 
malism developed earlier. Using equation (17), the local 
thermodynamic potential of the non-condensate atoms 
can be written as, 

1 h 2 
F(r) = -[g 5/2 (z eff ) ~ —V 2 V eff g 1/2 {z- ef f)], (27) 

where z eff = e^-^ff^ and V eff (r) = V(r) + 2up(r). 
For perturbative expansion we write the effective fugacity 
in the following form. 



z ef f = zexp[(3(ei - V(r) - 2up c (r) 

+2u( Pnc {0) - Pnc {r)))] (28) 

where, z — e^~ e °\ which is always less than 1. Assum- 
ing weak coupling, we expand the functions in terms of 
2u{p nc (Q) — p nc (r)). The perturbative expansion in this 
method turns out to be well behaved ||. Excluding the 
non-condensate part, the effective fugacity can be written 
as, 

z eff = ze/WrJ-ei) f or r<ro 

= ze ^~ v ^ for r > r , (29) 

where r$ is the turning point of Thomas-Fermi density 
p c . From the above expression we can see there is no ex- 
pansion parameter for condensate part. We approximate 
the local chemical potential within the condensate by its 
mean value. To calculate the free energy of the system 
within mean field theory, we subtract self energy from the 
total energy of the system. The self energy contribution 
is given by, 



E, 



u / P 2 nc (r)d 3 r. 



(30) 



In the following calculation we replace the coupling by 

a dimensionless parameter 77 = ( y/^jiV 1 / 6 ) 2 / 5 . Upto 

order ry 5 / 2 and Ti 2 the thermodynamic potential is given 
by, 

q(z,T)/N = t 3 g(4, z, 1/2, ar)0) + 
-^t 3 ^rj)[g 7/2 (z) + ~ar,0 g^- 2 ' 5 ^)] 

+2 V 5 / 2 t 7 / 2 [C(3/2)g 3 (S) - if(3/2, 3/2,3/2, 2)] 



f 
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t[g(2,z,l/2,a V p) + 



. - v W)b3/a(2) " -^hil2{ze- 2 l^)\ 

+2rf'H 1 ' 2 {am9i{i) ~ F(l/2, -1/2, 5/2, z))], (31) 
where the new functions and parameters are defined by, 



g(s,x,t,y) = —I(t,iy) 

i=l 

1 f°° 

I(t,y) = -= / dxe- x {x + ^ ( - t - 1 ') 
V 71 " Jo 



OO OO 



F{a,0, 1 ,x) = ^^- 



{x + y) 
x (i+j) 



a=\{l^N /N) 2 ^ 
t = 1/0 = kT/(hujN 1/3 ) 



(32) 

(33) 
(34) 

(35) 
(36) 



For small value of arj, the function <?(s, z, 1/2, ar]0) can 
be written as g s (z). Condensate fraction can be calcu- 
lated from the thermodynamic potential, 
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N ,N = l- Z 9 M^m (37 ) 
oz 

For T < T c , z w 1. To regulate some divergent terms 
appearing in the h 2 corrections, we set z w e P huJ , which 
is the natural energy gap from ground state. From the 
above expression we can see that the condensate frac- 
tion and other thermodynamic quantities are not fully 
scale invarient with respect to the dimensionless scaling 
parameters 77 and t. But the extra terms coming from 
the quantum corrections contain a factor 1/iV 2 / 3 . These 
terms depend on the number of particles in the system, 
and give finite size corrections to the scaling form of the 
thermodynamic quantities. In the large N limit and also 
in the high temperature phase, quantum corrections are 
negligible. But for small quantum confined systems they 
give finite corrections to the free energy. 




FIG. 1. Condensate Fraction for 10,000 atoms, 7/=. 18. 



the variation of condensate fraction with reduced tem- 
perature t, by choosing N — 10, 000 atoms and the di- 
mensionless parameter rj = 0.18. In Fig. 2 we have shown 
the finite size correction to the condensate fraction for 
same values of the parameters. From the graph we can 
see that the first order quantum correction is two order 
smaller than the leading order contribution. 

In conclusion, we have developed an extended Thomas- 
Fermi method at finite temperature for a system of 
Bosons. Systematic quantum corrections may be written 
in a series of a dimensionless parameter (3hu, which con- 
verges for the temperatures larger than the average spac- 
ing between energy levels. We have applied our method 
to a system of weakly interacting Bose gas confined in 
isotropic harmonic oscillator potential. Thermodynamic 
potential and condensate fraction are calculated analyti- 
cally upto if'/ 2 and H 2 order. Leading terms in H are the 
functions of scaling parameters t and 77 only, but leading 
order quantum corrections give N dependent corrections 
to the scaling form. In present situation quantum cor- 
rections are negligibily small. But quantum corrections 
will be important for small number of particles and for 
low dimensional systems. Also its magnitude depends on 
the nature of confinement potential. This method may 
be applied to investigate the properties of a system of 
confined charged bosons in weak magnetic field. 

I would like to thank M.V.N. Murthy and R.K. 
Bhaduri for helpful discussions. 




FIG. 2. Finite size correction to condensate fraction for 
10,000 atoms, where t is reduced temperature. 



For real physical systems the typical value of the pa- 
rameter 77 which determines the two-body interaction 
strength varies from 0.1 — 0.25. In Fig.l we have shown 
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